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Cross-checking the symmetry energy at high densities 
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By considering both the effects of the nucleon-nucleon short-range correlations and the isospin- 
dependent in-medium inelastic baryon-baryon scattering cross section in the transport model, two 
unrelated Au -|- Au experimental measurements at 400 MeV/nucleon beam energy are simultane¬ 
ously analyzed, a mildly soft symmetry energy (I/(po) = 37 MeV) at supra-saturation densities is 
obtained. This result is compatible with the recent result in Phys. Rev. C 92 , 064304 (2015) by 
comparing the available data on the electric dipole polarizability with the theoretical predictions. 
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I. INTRODUCTION 

The nuclear symmetry energy describes the single nu¬ 
cleonic energy of nuclei or nuclear matter changes as one 
replaces protons in a system with neutrons. Besides its 
impacts in nuclear physics [H, [2|, in a density range of 0.1 
^10 times nuclear saturation density, the symmetry en¬ 
ergy determines the birth of neutron stars and supernova 
neutrinos Q, a range of neutron star properties such as 
cooling rates, the thickness of the crust, the mass-radius 
relationship, and the moment of inertia The nu¬ 

clear symmetry energy also plays crucial role in the evo¬ 
lution of core-collapse simernova Q and astrophysical r- 
process nucleosynthesis Q. Thus the better we can con¬ 
strain the symmetry energy in laboratory measurements, 
the more we can learn from astro-observations. 

To constrain the symmetry energy in broad density re¬ 
gions, besides the studies in astrophysics [lol - [T^ . many 
terrestrial experiments are being carried out or planned 
using a wide variety of advanced new facilities, such as 
the Facility for Rare Isotope Beams (FRIB) in the US 
[T^ . or the Radioactive Isotope Beam Facility (RIBF) 
in Japan [l^ . To unscramble symmetry energy related 
experimental data, various isospin-dependent transport 
models are frequently used to probe the swnmetry en¬ 
ergy below and above saturation density |l|, 0. With 
great efforts, the nuclear symmetry energy and its slope 
around saturation density of nuclear matter from 28 anal¬ 
ysis of terrestrial nuclear laboratory experiments and as- 
trop hysical observations have been roughly pinned down 
[1^ , while recent interpretations of the FOPI and FOPT 
LAND experimental measurements by different groups 
made the symmetry energy at supra-saturation densities 
fall into chaos |16l - l^ . It does not seem to be clarified 
why the nuclear symmetry energy at supra-saturation 
densities is so uncertain, maybe the effects of pion in¬ 
medium effects [23l - l^ . the isospin dependence of in¬ 
medium nuclear strong interactions [2^, the short-range 
tensor force [13, are some factors. 

Recently, the high-momentum transfer measurements 
showed that nucleons in nucleus can form pairs with 
large relative momenta and small center-of-mass mo¬ 
menta [13, [30l| . This phenomenon was explained by the 
short-range nucleon-nucleon tensor interaction |3ll . l32j| . 


Such nucleon-nucleon short-range correlations (SRC) in 
nucleus lead to a high-momentum tail (HMT) in the 
single-nucleon momentum distribution [33l - l36l| . More in¬ 
terestingly, in the HMT of nucleon momentum distribu¬ 
tion, nucleon component is evidently isospin-dependent. 
The number of n-p SRC pairs is about 18 times that of 
the p-p and n-n SRC pairs [37l| . And in neutron-rich nu¬ 
cleus, proton has a greater probability than neutron to 
have momentum greater than the nuclear Fermi momen¬ 
tum [ 3 ^ . 

Unfortunately, effects of the above isospin-dependent 
SRC were seldom taken into account in most of currently 
used isospin-dependent transport models, while the lat- 
ters have been frequently used to unscramble symmetry 
energy related experimental data [ 13 , [ 13 , l33 - l43| . To ex¬ 
tract information on the symmetry energy from experi¬ 
mental data, in this study, by considering both the ef¬ 
fects of the isospin-dependent SRC and the important 
but often-overlooked in-medium baryon-baryon inelastic 
cross section in the isospin-dependent transport model, 
two unrelated experimental measurements are simultane¬ 
ously re-analyzed. 

II. THE IBUU TRANSPORT MODEL 

To probe the symmetry energy from experimental 
data, we use our recent updated Isospin-dependent 
Boltzmann-Uehling-Uhlenbeck (IBUU) transport model 
[i^l . In this IBUU model, nucleon-density distribution is 
given by 

r = i?(a;i)^^^; cosO = 1 — 212 ; (f = 2 ^ 3 ; 

X = rsinOcoscj)] y = rsinOsincj); z = rcosO. ( 1 ) 

Where R is the radius of nucleus, xi,X 2 ,x^ are three in¬ 
dependent random numbers. Since there is a depletion of 
nucleon distribution inside the Fermi sea, the proton and 
neutron momentum distributions with high-momentum 
tail reaching about 2 times the Fermi momentum [^ are 
given by the extended Brueckner-Hartree-Fock (BHF) 
approach by adopting the AV 18 two-body interaction 
plus a microscopic Three-Body-Force (TBF) [33|. Fig. [T] 
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FIG. 1: (Color online) Momentum distributions of neutron 
and proton in nucleus ^®^Au calculated with the BHF with 
Avl8+TBF. Taken from Ref. Q. 
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FIG. 2: (Color online) Evolution of the momentum distribu¬ 
tion of nucleon in nucleus ^®^Au. The upper panel (a) is for 
the case of initial nucleon distribution with HMT while the 
lower panel (b) is a simple Fermi-Dirac initial distribution. 
Both cases are under the interaction given by Eq. ©. 


shows such nucleon momentum distribution with high- 
momentum tail in ^®^Au. Compared with the distribu¬ 
tion in ideal Fermi gas, the excess energy of nucleon in 
colliding nuclei is subtracted from the total energy of re¬ 
action system. 

How long will survive the shape of the initial distribu¬ 
tions? To answer this question, Fig. [5] shows the plots 
of the implemented distributions at several times (0.5 
fm/c, 10 fm/c, 20 fm/c, 30 fm/c) for nuclei in the ground 
states when they are left to evolve under the interaction 
given by Eq. (© in this text and initialized in agreement 
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FIG. 3: (Color online) Time dependence of the number of 
evaporated nucleons by the nucleus in ground state using two 
initial distributions with and without high momentum tail as 
shown in Fig. [21 



FIG. 4: (Color online) Time dependence of the asymmetry 
Sbound = {pn - pp)l{pn + Pp) of bound nucleons in the nu¬ 
cleus in ground state using two initial distributions with and 
without high momentum tail as shown in Fig. [2] and Fig. (3] 


with the distribution seen in Fig. [TJ The lower panel (b) 
shown in Fig. [2] is under the simple Fermi-Dirac distri¬ 
bution as comparison. It is seen that the shape of the 
initial distribution of nucleon in momentum space is rel¬ 
atively well kept at the initial stage of collision at 400 
MeV/nucleon incident beam energy. Related to the pre¬ 
vious point, it is important to show a comparative study 
of the time dependence of the number of evaporated nu¬ 
cleons by the nuclei in ground state, both when are con¬ 
sidered a simple Fermi-Dirac distribution and a distribu¬ 
tion with high momentum tail. Shown in Fig. [H] is the 
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FIG. 5: (Color online) Kinetic symmetry energy and density- 
dependent symmetry energy with different x parameters. 
Note here that in heavy-ion collisions at 400 MeV/nucleon 
beam energy, the high-density part of the symmetry energy 
plays major role. 


time dependence of the number of evaporated nucleons 
by the nucleus in ground state using two initial distribu¬ 
tions with and without high momentum tail. One can 
clearly see that compared with the case without HMT, 
more nucleons emit from nucleus with HMT at the ini¬ 
tial stage of collision at 400 MeV/nucleon incident beam 
energy. However, to study the effect of the symmetry en¬ 
ergy, the asymmetry {Sbound = (p™ -Pp)/(Pn + Pp)) of the 
colliding system is more crucial. Shown in Fig. |T] is the 
time dependence of the asymmetry of bound nucleons in 
the nucleus in ground state using two initial distributions 
with and without high momentum tail as shown in Fig.[^ 
and Fig. [31 It is seen that with HMT the asymmetry of 
colliding nucleons in the reaction system is almost the 
same as that without HMT at the the initial stage of col¬ 
lision. One can thus conclude that although the present 
HMT considerations in initial colliding nuclei cause the 
instability of nuclei, it in fact does not affect much the 
present study of the symmetry energy at high densities 
using heavy-ion collisions. 

In this model, the isospin- and momentum-dependent 
mean-field single nucleon potential is used i.e., 


U{p,6,p,t) = Au{x)—+ Ai{x) — 

Po Po 

— - x6‘^) - 8xT^—^^—^6pr 

Po cr + 1 Po 


2Cr 


Po 


(fp 


fr{r,p') 


l + (p-p')2/A2 




po 


1 + {p-p'Y/hP‘ 


where r,r' = I/2(—1/2) for neutrons (protons), 5 = 
(pn - Pp)l{Pn + Pp) is the isospin asymmetry, and Pp 
denote neutron and proton densities, respectively. Specif¬ 
ically, the parameter values Au{x) = 33.037 - 125.34a: 
MeV, Ai{x) = -166.963 -b 125.34a: MeV, B = 141.96 MeV, 
Cr,r = 18.177 MeV, Cr,r' = -178.365 MeV a = 1.265, and 
A = 630.24 MeV/c are obtained by fitting seven empirical 
constraints of the saturation density po = 0.16 fm“^, the 
binding energy Eq = -16 MeV, the incompressibility Kq 
= 230 MeV, the isoscalar effective mass m* = 0.7m, the 
single-particle potential = 75 MeV at infinitely large 
nucleon momentum at saturation density in symmetric 
nuclear matter, the symmetry energy S{p) — 30 MeV and 
the symmetry potential = -100 MeV at infinitely 

large nucleon momentum at saturation density, frix^p) 
is the phase-space distribution function at coordinate r 
and momentum p and solved by using the test-particle 
method numerically (^ . Different symmetry energy’s 
stiffness parameter x can be used in the above single nu¬ 
cleon potential to mimic different forms of the symmetry 
energy. Since the kinetic symmetry energy, even its sign, 
is still controversial [46| , we at present give it a null value 
[ 4 ^ . For its density-dependence, we use similar form as 
that from the ideal Fermi gas model. Thus the density- 
dependent kinetic symmetry energy is expressed as 

= 12.5[(p/po)"/' - 1]. (3) 

Fig. [5] shows the kinetic symmetry energy we used and 
the density-dependent symmetry energy with different x 
parameters. It is seen that our used density-dependent 
kinetic symmetry energy is similar to that in Ref. (4^ . 

We can also see that x = 1,0,—1 cases roughly 
correspond positive slopes {L{pq) = 3pQdEsym{p)/dp) 
37, 87, 138 MeV, respectively. The following table shows 
the parameters used in Eg . ([2]) with SRC compared with 
the case without SRC [4^ . 


TABLE I: The parameters (in MeV) used in Eq. ([2]) with SRC 
compared with the case without SRC [4^ . 


jpkin 

-^sym 

0 ( = with SRC) 

12.5 (p/po)^^^ ( — without SRC) 

Aulx) 

33.037 - 125.34a; 

-95.98 - 91.157® 

Ai{x) 

-166.963 + 125.34a: 

-120.57 + 91.157® 

B 

141.96 

106.35 

a,r 

18.177 

-11.7 

Cr,r' 

-178.365 

-103.4 

a 

1.265 

1.333 

A 

630.24 

260 


The isospin-dependent baryon-baryon {BB) scatter¬ 
ing cross section (elastic or inelastic, including NN—i-NN, 
NA^NA, AA^AA, NN^NA) in medium cr^^dmm jg 

reduced compared with their free-space value by a 

factor of [ 4 ^ 


^medium (P; P) 


^medium /free 

B^elastic,inelastic ' ^^elastic,inelastic 

(p-bb/P'Bb)'^ , (4) 
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FIG. 6: (Color online) Reduced factor R^^dium of neutron- 
proton scattering cross section as a function of density and 
nucleonic momentum in asymmetric medium S = 0.2. The 
red shadow region and label ‘0.5~0.7’ denote experimental 
density and momentum region and the corresponding value 
of reduced factor 


where ^ibb and are the reduced masses of the col¬ 
liding baryon-pair in free space and medium (in medium, 
the effective mass of baryon is used), respectively. The 
effective mass of baryon in isospin asymmetric nuclear 
matter is expressed by 


m*g _ L tn-B dUs 
ms 1 p dp 


(5) 


For the baryon resonance A potential, the forms of 


= C/„, (6) 

+ (7) 

Ub = -^Un + -C7p, (8) 

= Up (9) 


are used [i^. As an example, the reduced factor 
^'medium neutron-proton scattering cross sec¬ 

tion in medium is shown in Fig. [B1 It is a function of 
density and nucleonic momentum in medium with asym¬ 
metry (5 = 0.2 for example. And generally R^edium > 
Rmed^um > K^edzum at Certain density and momentum. 
The nucleon-nucleon scattering cross section is reduced 
much in medium for colliding pair at high density and 
low momentum while it is less reduced at low density 
and high momentum. It is noticed from Fig. |6] that the 
used in-medium neutron-proton scattering cross section 
fits the experimental hard photon measurements quite 
well [^ . 


III. RESULTS AND DISCUSSIONS 



X 


FIG. 7: (Color online) Multiplicity of charged pion meson pro¬ 
duced in Au-|-Au reaction at 400 MeV/nucleon with different 
symmetry energies. The shadow region denotes the FOPI 
data |5l|. The dashed lines show the cases without SRG. 

Before studying the tt~ Jtt^ ratio, it is instructive to 
first see the production of charged pion meson in cen¬ 
tral An + An reaction at 400 MeV/nucleon beam en¬ 
ergy. Fig. [3 shows numbers of charged pion produced 
with different symmetry energies. It is seen that for x 
= I and X = 2 cases both produced tt~ and 7r+ fit the 
FOPI experimental data quite well. With stiffer symme¬ 
try energy a; = 0, the model gives somewhat smaller tt~ 
number than experimental data. While with very soft 
symmetry energy x = 3, the model gives both larger tt’*' 
and 7r“ numbers than experimental data. Comparing 
upper panel with lower panel, it is seen that sensitivity 
of the number of produced tt~ to the symmetry energy 
is at least 3 times that of 7r+. This is because the 7r“ 
mesons are mostly produced from neutron-neutron col¬ 
lisions, thus more sensitive to the isospin asymmetry of 
the reaction system and the symmetry energy [s^l . From 
this figure, one can see that without SRC, numbers of 
both 7r+ and tt~ are smaller than those with SRC. This 
is understandable since the SRC of nucleons increases the 
kinetic energy of nucleon. 

To reduce the systematic errors, most of the observ¬ 
ables proposed so far use differences or ratios of isospin 
multiplets of baryons, mirror nuclei and mesons, such as, 
the neutron/proton ratio of nucleon emissions, neutron- 
proton differential flow, 7r“/7r+, etc. Fig. [8] shows the 
7 r“/7r+ ratio predicted by our IBUU model with differ¬ 
ent symmetry energies. Because softer symmetry en¬ 
ergy causes more neutron-rich dense matter and 7r“’s are 
mainly from neutron-neutron collision whereas tt+’s are 
mainly from proton-proton collision , it is not surpris¬ 
ing that one sees larger 7r“ /tt'^ ratio with softer symme¬ 
try energy. To see the effects of SRC of nucleon-nucleon 
and the reduction of the in-medium inelastic baryon- 
baryon scattering cross section, with same x parame- 
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FIG. 8: (Color online) 7r“/7r^ ratio in Au+Au reaction at 400 
MeV/nucleon with different symmetry energies. Also shown 
are the effects of the SRC of nucleon-nucleon and the in¬ 
medium inelastic cross section on the ratio with same 

X parameters. 


FIG. 9: (Color online) The effects of in-medium inelastic cross 
section on the productions of different A resonance in Au-fAu 
reaction at 400 MeV/nucleon with the symmetry energy pa¬ 
rameter a; = 1. 


ters, we made calculations by turning off the SRC and 
the reduction of the in-medium inelastic baryon-baryon 
scattering cross section, respectively. From Fig. HI we 
can see that both of them affect the value of 7r“/7r+ 
ratio evidently. Both the SRC of nucleon-nucleon and 
the reduction of the in-medium inelastic baryon-baryon 
scattering cross section decrease the value of 7r“/7r+ ra¬ 
tio evidently. Proton-proton collision is also affected by 
the Coulomb action, so tt'*' production, which is mainly 
from proton-proton collision, is relatively less affected by 
the reduction of the in-medium inelastic baryon-baryon 
scattering cross section. Fig. shows the evolution of A 
resonance production with free and in-medium inelastic 
cross sections. It is clear seen that the in-medium inelas¬ 
tic cross section affects the production of A“ (decay into 
TT~) much than A++ (decay into 7r+). This is the reason 
why the reduced in-medium inelastic baryon-baryon cross 
section decreases the value of 7r“/7r+ ratio. The SRC of 
neutron and proton causes small asymmetry of matter, 
which corresponds small value of 7r“/7r+ ratio. Therefore 
both the SRC of nucleon-nucleon and the in-medium in¬ 
elastic baryon-baryon cross sextion should be taken into 
account in transport calculations. From Fig. HI we can see 
that the FOPI pion experimental data supports a softer 
symmetry energy (x = 1, 2, even x = 3). Note here that 
the density region probed here is about 1 ~ 1.5 times 
saturation density [^. So we do not conclude what are 
the constraints of the values of the symmetry energy and 
its slope around saturation density. 

While A and tt’s scattering and re-absorption can de¬ 
stroy the high density signal in a certain degree. Treat¬ 
ment of Delta dynamics in transport models is not so 
straightforward such as the competing effects of the mean 


fields and A thresholds. To understand quantitatively 
the symmetry energy effect on pion production, it is im¬ 
portant to include the isospin-dependent pion in-medium 
effects [25 - E^ . And recent work of MSU group 
demonstrates that the ratio of tt’s spectra is more sensi¬ 
tive than the ratio of integrated yields because the latter 
gives ambiguous result since it does not distinguish tt’s 
messenger of high density from the rest. Therefore, more 
theoretical and experimental studies are needed to pin 
down the high-density behavior of the symmetry energy 
by pion probe. 
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FIG. 10: (Color online) Elliptic flow of emitting nucleons in 
Au+Au collision at 400 MeV per nucleon incident beam en¬ 
ergy with different symmetry energies. The shadow region 
denotes the experimental FOPI-LAND data [ 2 ^ . The dashed 
lines show the cases without SRG. 

To cross-check the symmetry energy over a broad re¬ 
gion (x = 1, 2, 3), one has to search for other constraints. 
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Fig. [To] shows predicted neutron and proton elliptic flows 
in An + An reaction under the FOPI-LAND experimen¬ 
tal conditions (^re we use the data with b = 7.2 fm case) 
and geometry [20| . The experimental data was multiplied 
by a factor 1.15 owing to dispersion of the reaction plane 
[^ . From Fig.[Tni(a), it is seen that our model give some¬ 
what lower value of the elliptic flow of neutrons. How¬ 
ever, from Fig.[Tni(b), it is seen that the predicted elliptic 
flow of protons fit experimental data quite well. In our 
model, free nucleons are identified by their local densities 
Pfree < Po/8, which Corresponding to deuteron’s nucleon 
average density O.Q2fm~^. The identification standard 
of free nucleons affects the value of V.p and V 2 , but does 
not affect the ratio of Vp/Vf much. With stiffer symme¬ 
try energy, the value of nucleon elliptic flow should be¬ 
come larger. However, this trend seems not right for the 
stiffer symmetry energy a; = 0. This abnormal behavior 
may be caused by the competing effects of the SRC and 
the symmetry energy and deserve further study. From 
Fig. 1101 at such experimental conditions and geometry, 
effects of the symmetry energy on both proton and neu¬ 
tron elliptic flows can not be seen clearly. One way to 
enlarge the effects of the symmetry energy is the relative 
changes of proton and neutron elliptic flows, such as the 
ratio of neutron and proton elliptic flows ■ From 

this figure, one can see that without SRC, values of both 
neutron elliptic flow and proton elliptic flow are larger 
than those with SRC. This is understandable since the 
SRC of nucleons decreases the anisotropic emissions of 
neutrons and protons. 



FIG. 11: (Color online) Same as Fig. 1101 but for the ratio 
of IV 2 ■ The effects of the SRC of nucleon-nucleon on the 
V 2 IV 2 with same x parameters are also shown. 

Shown in Fig.fTTlis predicted elliptic flow ratios of neu¬ 
tron and proton V 2 with different symmetry energies 
as well as experimental data . Since stiffer symmetry 
energy/symmetry potential causes relatively more neu¬ 
trons to emit in the direction perpendicular to the re¬ 


action plane [H^, one sees larger values of elliptic flow 
ratios of neutron and proton V 2 IV 2 with stiffer symme¬ 
try energies. With the SRC of nucleon-nucleon in the 
transport model, values of the 1^2 ratio are larger 
than that without the SRC of nucleon-nucleon. This is 
because the SRC of nucleon-nucleon cause neutron and 
proton to be correlated together, the value of 1^2 ra¬ 
tio trends to unity. Owing to the competing effects of the 
SRC and the symmetry energy, for a: = 0 case, the effects 
of symmetry energy on the trend of the ratio of 
with the SRC changes compared with that without the 
SRC. 

On the whole, the sensitivity of the observable V 2 IV 2 
to the symmetry enegy at FOPI-LAND experimental 
conditions and geometry is smaller than that of the FOPI 
7 r“/7rT ratio. Other nucleon observables should be fur¬ 
ther explored. 

Fig. Ill] indicates the FOPI-LAND elliptic flow experi¬ 
mental data does not favor very soft symmetry energy {x 
= 2, 3). Combining the studies of nucleon elliptic flow 
and previous tt~ ratio, one may roughly obtain the 
symmetry energy stiffness parameter a; = 1. It in fact 
corresponds a mildly soft density-dependent symmetry 
energy at supra-saturation densities. While the specific 
density region of the present constraints on the nuclear 
symmetry energy needs to be further studied [s^. 

The small effects of the symmetry energy on pionic 
and nucleonic observables, pion suffering from unclear 
TT — N — A dynamics and pion in-medium effect, neutron 
detection efficiency, bound or unbound nucleon identifi¬ 
cations, nucleon in-medium and isospin strong interac¬ 
tions as well as all kinds of experimental measurement 
errors, etc., all affect the probe of the symmetry energy, 
not to mention uncertainties and complexities of nuclear 
transport models, thus fully convincing constraints of the 
symmetry energy at high-density are not easy to achieve. 

It is, however, interesting to see that the present result 
on the symmetry energy stiffness parameter a; = I (which 
corresponds a slope of symmetry energy at saturation 
density (L(po) = ipodEsym{p)/dp = 37 MeV) agrees 
with the recent result L{po) = 20 - 66 MeV quite well 
by comparing the available data on the electric dipole 
polarizability with the predictions of the random-phase 
approximation, using a representative set of nuclear en¬ 
ergy density functionals (^ . 


IV. CONCLUSIONS 

In summary, by incorporating the short-range corre¬ 
lations of nucleon-nucleon and the in-medium inelastic 
baryon-baryon scattering cross section into the isospin- 
dependent transport model and based on the FOPI 
and FOPI-LAND experimental measurements, I cross¬ 
checked the ratio and the ratio of neutron ellip¬ 

tic flow and proton elliptic flow Vp/V 2 in Au-I-Au colli¬ 
sion. A mildly soft symmetry energy at supra-saturation 
densities supports both FOPI and FOPI-LAND experi- 
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mental measurements. The studies also show that both 
the short-range correlations of nucleon-nucleon and the 
in-medium inelastic baryon-baryon cross section play im¬ 
portant role in probing the symmetry energy with heavy- 
ion collisions. 

Since the symmetry energy plays crucial roles in both 
nuclear physics and astrophysics, more subjects, such 
as the density region that some observables probed and 
more sensitive observables to the symmetry energy at 
high densities, deserve further study. 
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